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Abstract 

The linear instability of the hypersonic boundary layer on a curved wall is considered. 
As a starting point the viscosity of the fluid is taken to be a linear function of temperature 
and real-gas effects are ignored. It is shown that the flow is susceptible to Gortler vor- 
tices and that they are trapped in the logarithmically thin adjustment layer in which the 
temperature of the basic flow changes rapidly to its free stream value. The vortices decay 
exponentially in both directions away from this layer and are most unstable when their 
wavelength is comparable with the depth of the adjustment layer. The non-uniqueness 
of the neutral stability curve associated with incompressible Gortler vortices is shown to 
disappear at high Mach numbers if the appropriate ‘fast’ streamwise dependence of the 
instability is built into the disturbance flow structure. It is shown that in the hypersonic 
limit wall-cooling has a negligible effect on the stability of a fluid with a given value of the 
Chapman constant. 


* This research was supported in part by the National Aeronautics and Space Admin- 
istration under NASA Contract No. NASl-18605 while the first author was in residence 
at the Institute for Computer Applications in Science and Engineering (ICASE), NASA 
Langley Research Center, Hampton, VA 23665. 
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1. Introduction 


Our concern is with the linear development of Gortler vortices in hypersonic boundary 
layers. Interest in the stability properties of hypersonic boundary layers has been stimu- 
lated by recent research and development into hypersonic aircraft. In some situations it 
is envisaged that Mach numbers as high as 20-25 will be achieved; in that situation real 
gas effects will certainly be important. Here, as a starting point, we ignore such effects 
and assume that the fluid under investigation behaves as a perfect gas. The other major 
simplification which we make is that the viscosity of the fluid is given by the Chapman’s 
viscosity law. This simplification can of course be justified but at any rate it is probable 
that our analysis could be extended to deal with more realistic viscosity-temperature laws. 

We shall take into account the effect of boundary layer growth on the vortices using 
asymptotic methods related to those used for the incompressible case by Hall (1982a, b, 
1983, 1985, 1988) and Hall and Lakin (1988). The main feature of the incompressible 
case is that the growth of vortices with wavelengths comparable with the boundary layer 
thickness is dominated by non-parallel effects and can only be described self-consistently 
by integrating numerically the linear partial differential equations obtained by perturbing 
appropriately the Navier-Stokes equations. These equations were solved by Hall (1983) 
who showed that the position of neutral stability of a Gortler vortex is a function of its 
upstream origin; this means that in this wavelength regime the concept of a unique neutral 
curve is not tenable. Thus the results of, for example, Gortler (1940), Hammerlin (1956) 
which were obtained by ignoring the streamwise dependence of the vortices are not valid. 
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However, at smaller vortex wavelengths the vortices only feel the local structure of 
the boundary layer and localize themselves so as to maximize their downstream growth. 
In the neutral case this position corresponds to the location where Rayleigh’s criterion is 
most violated. Moreover the vortices then develop downstream in a quasi-parallel manner 
and Hall (1982a) showed how the Gortler number could then be simply expressed as an 
asymptotic expansion in terms of the wavenumber. An asymptotic investigation of the 
equation used by the parallel flow theories shows that they give the first few terms in 
this expansion but that they are in error when non-parallel effects become important. 
Thus, in the only regime where the parallel flow theories are valid, a relatively trivial 
asymptotic expansion of the Gortler number is available and is indeed more accurate than 
that obtained by a numerical solution of the parallel flow eigenvalue problem. 

In a recent paper Hall and Malik (1987) applied the approach of Hall (1982a) to 
compressible boundary layers at O(l) Mach numbers. A parallel flow theory of Gortler 
vortices in compressible boundary layers at 0(1) Mach numbers had previously been given 
by El-Hady and Verma (1981) who followed the approach of Floryan and Saric (1979). 
Here, as a first step towards an understanding of the structure of Gortler vortices at 
hypersonic speeds, we look at the limiting form of the Hall-Malik calculation at high Mach 
numbers. 

It turns out that the “most dangerous” location in a hypersonic boundary layer is 
the logarithmically small region where the basic state temperature field adjusts to its free 
stream value. It is in this layer that asymptotically small wavelength vortices described 
by the structure found for the incompressible case by Hall (1982a) become trapped. 
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The main result of the present calculation is that the logarithmically small layer at 
the edge of a hypersonic boundary layer can support a range of Gortler vortex wavelengths 
other than those appropriate to the high wavenumber limit of Hall (1982a). Moreover the 
most dangerous wavelengths of Gortler vortices at hypersonic speeds are comparable with 
the thickness of the logarithmic layer. Furthermore, the streamwise development of the 
vortices can be taken care of by using a similarity variable and a multiple scale or WKB 
approach so that the structure of the vortices is obtained by solving ordinary differential 
equations. However, it should be stressed that the replacement of the streamwise derivative 
of the disturbance by a constant in the work of El-Hady and Verma means that their results 
are in error. Thus at hypersonic speeds Gortler vortices have a quasi-parallel behaviour 
and a structure considerably simpler than that relevant to the 0(1) Mach number regime. 
We shall see that the neutral curve at hypersonic speeds is well-defined with left and right 
hand branches and a unique minimum Gortler number above which linear perturbations 
will grow. 

The procedure adopted in the rest of the paper is as follows: in §2 we formulate 
the linear stability problem for a compressible boundary layer. We then write down the 
eigenrelation found by Hall and Malik (1987) and by taking the further limit of large 
Mach number show the importance of the logarithmically thin layer at the edge of the 
boundary layer. In §3 we show that this layer can in fact support a whole range of vortex 
wavenumbers and show that there exists a unique critical Gortler number at hypersonic 
speeds. In §4 We present results for a variety of flow properties and draw some conclusions. 
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2. Formulation of the stability equations and the small wavelength hypersonic 
limit. 

Our formulation follows closely the derivation of El-Hady and Verma (1981) who 
extended the usual incompressible approach in a straightforward manner. Suppose that 
L is a typical streamwise length scale while U 0 Q , p^, T^, Hqq are the free stream values 
for velocity, density, temperature and viscosity. If we take ( x*,y*,z *) to be a coordinate 
system with x* measuring distance along a rigid wall of variable curvature -j/c( ^-) then 
we define a curvature parameter 6 by 

f 2 - 1 ' 

and consider the limit 6 — > 0 with the Reynolds number R defined by 

R = C/ ° ojLPo ° , (2.2) 

Moo 

taken to be so large that the Gortler number 

G = 2R*6 (2.3) 


is 0(1). We define a Mach number M by 


Af 2 = 


uL 

iR*T 0 o ■ 


(2.4) 


Here R* is the gas constant whilst 7 is the ratio of the specific heats of the gas. The 
coordinate y* measures distance normal to the body whilst z* measures distance in the 
spanwise direction .The reader is referred to the book by Stewartson(l964) for a more 
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detailed discussion of compressible boundary layers. We assume that the basic flow corre- 
sponds to a zero pressure gradient and that the viscosity /x* and temperature are related 
by Chapman’s Law: 

f -of. 

Moo ■*- oo 

Here C is a constant and T* is the temperature. We denote the Prandtl number by T, 
and define dimensionless variables x, y, z by 

(x, y, z) = {x*L-\y*R±L-\z*R l *L- 1 ), 

whilst (u,v,tv), scaled on U Q Q , R~^Uoo, and R~*Uoo, denote the velocity corresponding 
to (x, y,z). Finally we let p, T and n denote density, temperature and viscosity scaled on 
their free-stream values. The basic state is then given by 

(u, r, id) = (u(x, y),v(z,y), 0), T = T(x,y), p = p(x,y), p = p (x,y). (2.6) 


( 2 . 5 ) 


If we define the Howarth-Dorodnitsyn variable y by 



then u and v can be written as 


u = f'{r)),v= -^=(r(r?)(7?/'-/) - f [ rjT'{ri)dr)), 
V2x Jo 


where 


V = 


y 

y/2x 


and 

/(„) = ^CF(J=), 


( 2 . 7 ) 


(2.8a) 


( 2 . 86 ) 
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where F is the Blasius function. The temperature field T(rj) is then written in the form 

f°° f°° /* r?1 

T = l + B {f"{vi)} r din + T{'i-l)M 2 {f"{rh)} r dT}i {/"(?72)} 2_r ^^2 - (2.9) 
Jff Jff Jo 

Here B is determined by the appropriate boundary condition to be imposed on T at the 
wall. In the aerodynamic situation it is expected that the wall will be cooled since otherwise 
the wall temperature will become intolerably large at high Mach numbers. However, in 
order to indicate the appropriate vortex structure at high Mach numbers we shall, in the 
first instance, assume that the Prandtl number is 1 and that the wall is adiabatic. In that 
case T reduces to 

T = 1 + - f' 2 ). (2.10) 

Next suppose ,as in Hall and Malik(1987), that the flow is perturbed to a spanwise periodic 
stationary vortex structure with wavenumber a. The linearized stability equations for 
Gortler vortices are then found by linearizing the Navier-Stokes equations about (2.7), 
(2.9) and retaining the leading order terms in the high Reynolds number limit. We obtain 

^{uU) x -^-p,a U-\-^VUy—(fiUy)y-\-^UyV—[-^(UU x -j-VUy)-i-(flUy)y}T—j/,UyTy = 0, (2.11a) 

^(u x + KuG)U- Cp.yU x — (c + 1) fiUxy — p, x U y + -^(vV)y + ^?V X + p,a 2 V - (c + 2) (/xVy) y + Py 

-[J^(uv x + VVy + ±kGu 2 ) + (c + l'jflUxy + CjlyU x + (c + 2) (j WVy)y + frzUy]T - u y T x 

— \cjxu x + (c + 2)jj.Vy\T y — cfiyiaW — (c + l)iap,W y = 0, (2.116) 
fi x iaU + (c + l)p,iaU x + p, y iaV + (c + l)p.iaV y — iaP + cfr(u x + v y ) 
iaT - ^ W x - (c + 2)fia 2 W - ^W y + (fi W y ) y = 0, (2.11c) 
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(y). + (y), + «(?) - (p), - (p)„ - (2. nd) 

fT,u - 2(7 - i)M*im v u„ + ir„K - [i(of, + «r„) + (7 - + i 

(AT„)„]r + isr, + ^ a2 T + (is - pf„)r # - i(pr„)» = o.( 2 .iie) 

Here /2 = dfi/dT and c = (A /ft) - 2/3 where A is the bulk viscosity, whilst ( U,V,W ) , 
P, and T denote the vortex velocity field , pressure and temperature respectively. In the 
calculations reported later in this paper c has the numerical value —.666666 

It was shown by Hall (1982a) that in the incompressible case the neutral curve for 
small wavelength vortices has G ~ a 4 and the vortices are confined to a layer of depth a~ s 
where the flow is locally most unstable. Hall and Malik (1987) extended this approach to 
the above system for M = 0(1) and wrote 

G = g 0 a 4 + (^a 3 H . (2.12) 


In the neutral case the vortices are trapped at the location (in y ) where 

£ * _ £ 2 , gou 2 ^ uu y g 0 _ A 

“ 2f r y T 2 r~ ’ 


(2.13) 


and 


d6* d 2 S n 

W = 0 ’ W <0 - 


(2.14) 


In fact (2.13) determines the constant go for a given value of the streamwise variable x 
such that the vortices are locally neutrally stable there, Wadey (1989) has discussed the 
strongly nonlinear compressible Gortler vortex regime using the approach of Hall and Lakin 
(1988) and shows that (2.13) then determines u such that a strongly nonlinear vortex can 
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exist. For simplicity we take T = 1, in which case (2.13) reduces, after some simplification, 


to 

JO cr»{i + to=?«L(i - /' 2 )> e 

/>/»{ i + ’ 1 ' ‘ 

where (2.5) has been used and a prime denotes differentiation with respect to 77 . In 

Figure 1 we have shown go/yj2x calculated from (2.15), subject to g' 0 = 0 , g" > 0 , as a 
function of Mach number with C = 1 , 7 = 1.4. We see that go initially increases from its 
incompressible value, reaches a maximum at Mo ~ 5, and then decays. The corresponding 
value of 77 where go satisfies (2.15) is denoted by rj c and is shown as a function of M in 
Figure 2. We see that 77 c moves into the free stream as M increases. The precise large M 
asymptotic structure of go and 77 c is found by first noting that for 77 >> 1 


f ~ y/d ^~ 0 + ( j / Vc -/?) 2 + "' } ’ 


where 0 = 1.2168 and c is a constant. It then follows that the minimum value of go for 
M » 1 occurs for 


2 

(n/Vc-fi) 


M~ 2 . 


More precisely if we define y by 


JL 

Vc 


— 0 = V 2 log M 2 — 


{y + log \Z2C log M 2 } 
^/2 log M 2 


we find that (2.15) reduces to 

go _ 2C 5 / 2 {1 + JVe^} 6 

\/2x Ne^V 2 logM 2 


(2.16) 
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with N = (7 — 1 )c\/C, so that go attains its minimum value when Ne& = ^ and then 

» - 10Cl (f) 8 (217) 

V2i V21 ogM*' y ' ’ 

and this asymptotic prediction is shown in Figure lby the dotted curve. In fact in the 

layer where y = 0(1) the temperature is given by 

T = 1 + Ne», (2.18) 


so that the vortices are trapped in the logarithmically thin layer where T adjusts from 
its free stream value of 1 to an 0(M 2 ) temperature typical of a hypersonic boundary 
layer. Thus in the double limit a — ► 00 , M — ► 00 the vortices are trapped in a thin layer 
embedded within the logarithmic layer. We now show in the next section that other more 
unstable vortices persisting throughout the adjustment layer can exist. 

3. Hypersonic Gortler vortices with wavelength 0{(\/21og M 2 ) X ) 

We now investigate the possibility that the logarithmic adjustment layer at the edge 
of the boundary layer can support vortices with larger wavelengths than those described 
by Hall and Malik. If the nondimensional wall temperature of the basic state is T w then 
B in (2.9) is determined by 


„ r„ - 1 - rh - 1 i)) r *>i Jo ” 1 (/"M ) 2 - r *! 2 

& ^OO ] 


(3.1). 


/o (/"(<Il) r *»l 

If we impose an adiabatic wall temperature condition then (3.1) determines the wall re- 
covery temperature T w and B — 0. Since the Prandtl number is not necessarily unity now 
we now define a stretched variable y in the adjustment zone by 


— Vcfl = y log A/ 2 j=L= . 


(3.2). 
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After some manipulation we can then show from (2.9) that for y = 0(1) the temperature 
T may be expressed as 

f = 1 + Ne» (3.3) 

where 

N^j(^=) r {^-l)T / o °V(»)) 2 - r A; + UmJBM-*)}. (3.4) 

Thus the wall condition on T determines the constant N appearing in (3.3) which gives 
the temperature in the logarithmic adjustment zone. In fact we shall see that the critical 
Gortler number at which instability first occurs is independent of N, hence in a hypersonic 
boundary layer wall cooling has no effect on the onset of the Gortler instability. 

For convenience we now define Ml and M2 by 

Mi = M 2 , M 2 = \/2 log M 2 , (3.5a, b) 


and if we are to seek a vortex structure with wavelength comparable with the depth of the 
adjustment layer we must expand the wavenumber a in the form 


M 2 k 

a = — + 


y/2x 


(3.6) 


Here the factor \/2x has been inverted in order to scale x out of the eigenvalue problem 
G = G(x , a); we are thus determining G such that the flow is locally neutral at the position 
x. The growth or decay rates of a vortex trapped in the adjustment layer must be 0(M 2 ) in 
order to enter the zeroth order stability problem so that in the hypersonic limit the parallel 
flow approximation becomes valid. Here we concentrate on finding the neutral curve and 
therefore determine G = G{x,a) such that the vortex with wavenumber a is neutrally 
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stable at x. The previous discussion of the double limit a — ► oo , M — ► oo suggests that 
the Gortler number G should be expanded as 


G = 


MiG 

(2x)f + 


However an examination of (2.11b) shows that the term -~^\kGu 2 T which balances with 
the term ua 2 V in the limit G — » oo will be comparable with the term ^2 uv x T when 
G = 0 (M 2 ). Thus , as pointed out by F.T. Smith (private communication) , the curvature 
of the basic state produces an effective negative Gortler number of order M 2 in the absence 
of curvature so that instability cannot occur for G « M 2 . Hence we must expand G in 
such a way that the dominant contributions from the terms uv x and £ GKu 2 cancel in the 

adjustment layer. To achieve this we note that in this layer uv x = -^ QM 2 H where 

Q = jtfT / 0 °° r\T’dy. Thus G must be expanded in the form 


G = 


1 

(2x) 3 / 2 


{QM 2 + MiG + •••}, 


(3.7) 


and the disturbance quantities expand as 


| 

U= [U{x,y) + ---]E,V = [M 1 M 2 V (x, y) -\ ]E,W = \M x M 2 W{x, y) + •••]£, (3.8a - e) 


P = [MiM 2 2 + • • •] E., T = [M x + • • •] E 

where E = e-f M2<T ( z ) dz . Thus the local growth rate of the vortex is M 2 cr(x) but here 
we restrict our attention to the situation a(x) = 0. Again the factors of \/2x inserted in 
(3.8) have been introduced in order to scale x out of the eigenvalue problem at the neutral 

! 

! 
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location. If the above expansions are inserted into (2.11) we find after some manipulation 
that 

fT dV de Net . IT- ~ - 
\lcij + dj = ~r^c v+ ^ + ,km ' 

[ TdP . . d 2 # << n ,TNey . [cdo 

\lcls ={c + 2) ' m W + {{c + 2) —- 1} yrTy 


dy 2 ' lv ~ ' 7 

VrcN 2 e 2 t. 


dW f Fc _ 

dy + Wr 7 


(c + 4) 


- fa* + (c + 2) — 1 j r— - ~ } - ,*/VCT 

\/!^JVe^}*W+{C7 2 P+^^-(c+2)r^}F,r^ = ^+CJ 2 P^(3.9a,6,c,d,e) 
+(c + 2) v^iVe W + tijP - (c + l)CJtfc^ + C Net ikd, ^ ^ J" 1 + { - iV 

dy ntt^ nit k P 


dy 2 dy * 


5 *(1 + 7- 1 )7~ 1 }d - ^yjLyJ = V = W = 0 , y = ±oo. 


Here the function 7{y) = 1 + Net. 

The zeroth order approximation to the x momentum equation has not been written 
down here because it decouples from the other equations and thus plays no role in the 
eigenvalue problem. A result of some importance follows immediately from (3.9) if we 
note that by a simple change of origin in y we can, without any loss of generality, set 
N = 1 since whenever it appears it is multiplied by et . Again this means that the wall 
conditions on the basic temperature field have no effect on the Gortler number for the 
onset of Gortler vortices in the hypersonic limit. 

The eigenvalue problem specified by (3.9) was solved using the compact finite differ- 
ence scheme of Malik, Chuang and Hussaini (1982). In order to apply the scheme the 
boundary conditions (3.9e) must be applied at finite values of y so it is necessary to in- 
vestigate the required decay of the disturbance quantities at y — ±oo in order to derive 
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asymptotic boundary conditions. After some manipulation it is found from (3.9) that when 


y — ► -oc 


0 = Pe^ 1 * V = Qe M3 *' + Re ^ H ^ 2 ^ — — — 111 — (3.10) 


where 


1, / 4Cfc2 11 /i 4Cfc 2 . r N 

mi = 2 {* v 1 ^ — p — } > ^2 — 2 + V r 2 ^ "* f — — p"‘( 3,lla ’ c ) 


The constants P, Q and R are arbitrary so that as y -* oo we have three independent 


solutions of (3.9). The quantities W, P, can then be expressed in terms of P, Q and 

R and when the three latter constants are eliminated we obtain the asymptotic boundary 
conditions 


tjt - Ml 9 = 0 ,ikW + F0 = t/J/x 2 V + {-Sfi 2 j^+ Eitl—tAye + ^1(M3 - M 2 )P, 
ay V o V o M 3 M 3 

+ = ^ + + + ^ f 


where 


S = = c (c ^ M r- 1 ). 

When y — ♦ 00 a related set of boundary conditions can be derived using the WKB 
method. We do not give the full details here because it is a routine but tedious calculation 
to derive these boundary conditions. However in order to point out the rather abrupt 
decay of the disturbance which occurs for y -* 00 we note that in this limit (3.9d) reduces 


%■ - 
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and if we set V =0 the required decaying solution of this equation has 

0 ~ e~V^7r*e s . 

In fact, if we do not set V =0 and consider the limit of the other disturbance equations we 
find two further decaying solutions with a similar structure to that above. The appropriate 
boundary conditions, as y —* +oo are then found by eliminating the constants associated 
with the three decaying solutions. However the structure of the decaying solutions for 
y — * oo means that the lower boundary of the region of vortex activity will be much 
sharper than the upper one. This difference will clearly be most pronounced in the small 
wavenumber limit. 

We shall present the results of our numerical investigation in the next section. Here we 
note that in order to find the neutral Gortler number correct to the graphical accuracy of 
the figures shown in that section it was sufficient to use 200 equally spaced grid points and 
approximate (— oo,+oo) by (—10., 15.0). Without the asymptotic boundary conditions at 
±oo the magnitude of the required approximation to — oo must, particularly for k « 1, 
be significantly increased. 

4. Results and Discussion 

We shall in the first place discuss the results of our numerical solution of the eigenvalue 
problem G — G[k). This was carried out using the scheme discussed in the previous section 
and Figure 3 shows the neutral curves for three different values of the Chapman constant C 
and T = 0.72. We see that the curves have well defined left and right hand branches typical 
of Taylor vortex problems. The critical Gortler number is seen to increase monotonically 
with the Chapman constant C. We note that for a given value of C the basic hypersonic 
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boundary layer is unstable to Gortler vortex disturbances if (k, G) is above the appropriate 
neutral curve of Figure 3. 

Thus we see that in the hypersonic limit a fluid with viscosity given by Chapman’s 
Law supports Gortler vortices whose growth is governed by a parallel flow theory. The 
parallel flow approach is valid in this situation because of the development of the thin 
logarithmic adjustment layer in a hypersonic boundary layer. At 0(1) Mach numbers 
there is no such adjustment region and the growth of the Gortler vortex mechanism is 
controlled by non-parallel effects. 

The eigenfunctions associated with the C = 0.5 calculations are shown in Figures 4 
and 5 for two values of the vortex wavenumber. We note that the rate of decay of the 
vortex when y — ► oo is much greater than when y — ► — oo. This is consistent with the 
discussion of the previous section. 

Now let us turn to a discussion of the relationship of our results with those found in 
earlier investigations of Gortler vortices in compressible boundary layers. Firstly, we note 
that the scalings used in the present paper are those inferred by the small wavelength limit 
of Hall and Malik (1987) if the further limit M -* oo is taken. However, the fact that 
the logarithmic adjustment layer is relatively thin means that both the left and right hand 
branches of the neutral curve can be described by a quasi-parallel theory. In actual fact 
it is perhaps strictly not correct to describe the analysis of §3 as a parallel flow theory 
because the whole basic flow and disturbance structure in the adjustment layer depends 
on x. In a parallel flow theory, e.g. El Hady and Verma (1981), the operator ud x + vd y , 
which appears in (2.11), would be replaced by f3u + vd y where /? is a constant growth 
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rate. In our analysis this operator is rewritten in terms of d x and dy with the result that 


the x derivative in the resulting operator becomes formally O { ( \/2 log M 2 ) ~ 1 } smaller 
than the y derivative and therefore negligible. In other words the disturbance has a ‘fast’ 
x dependence built in when we change variables; this structure is lost when a parallel 
flow theory of the type used by El Hady and Verma arbitrarily replaces ^ acting on a 
disturbance quantity in (2.11) by a constant. 

El Hady and Verma gave neutral curves for Mach numbers in the range 0 < M < 5 
and found that the critical Gortler number increased monotonically with Af. (See Figure 
1 of El Hady and Verma (1981)). This is entirely consistent with our prediction that the 
Gortler number scales on (log M 2 ) * for M » 1. However in our calculations we find that 
the critical wavenumber increases like (log M 2 )* whilst they found that , at least in the 
range 0 < M < 5, it decreased with M. 

The first possible explanation for this difference is that El Hady and Verna did not 
perform computations at large enough values of M in order to see the correct large M 
dependence of the wavenumber. Alternatively, it could be that the parallel flow approx- 
imation in their calculation necessarily leads to the incorrect wavenumber behaviour for 
large M. At 0(1) values of M their results are clearly incorrect and the disturbance growth 
must be described using the approach of Hall (1983). In any case it is clear that in the 
only Mach number regions where a parallel flow theory is valid the approach of El Hady 
and Verna effectively ignores the ‘fast’ x dependence of the disturbance in the adjustment 
layer.The authors acknowledge the comments made by the referees on the original form of 
this paper. 
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Legends 


Figure 1. The dependence of g 0 /y/2x on M calculated from (2.15) 

Figure 2. The dependence of r) c , the value of r\ at which go is calculated, on M 
Figure 3. The neutral curves in the G — k plane for three different values of the 
Chapman constant C. 

Figure 4. The eigenfunctions for the case C = .5,fc = .4 
Figure 5. The eigenfunctions for the case C = .5, k = .2 
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